Proposal for a Nuclear Gamma-Ray Laser of Optical Range. 
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A possibility of amplification of the 7.6 eV 7-radiation by the stimulated 7-emission of the en- 
semble of the 229 ™Th isomeric nuclei in a host dielectric crystal is proved theoretically. This ampli- 
fication is a result of: 1) the excitation of a large number of the 229m Th isomers by laser radiation; 
2) the creation of the inverse population of nuclear levels in a cooled sample owing to the inter- 
action of thorium nuclei with the crystal electric field or with an external magnetic field; 3) the 
emissions/absorption of the optical photons by thorium nuclei in the crystal without recoil; 4) the 
nuclear spin relaxation through the conduction electrons of the metallic covering. 

PACS numbers: 23.20.Lv, 27.90.+b, 42.72.Bj 



Nuclear 7-ray laser is a known challenge of the mod- 
ern physics. To succeed in the field, one has to overcome 
at least two basic problems, namely: to accumulate a 
sufficient amount of isomeric nuclei in a sample, and to 
narrow down the emission 7-ray line to its natural radia- 
tive width. These and other problems, as well as early 
unsuccessful attempts of their solution are discussed in 
detail in reviews (see also list of reviews in 

Nevertheless, there is a plausible way to overcome all 
difficulties and to develop a unique 7-ray laser working on 
magnetic dipole (Ml) transition in the optical range be- 
tween the first excited level (Jf s = 3/2+) and the ground 
state (JJ r = 5/2+) of the 229 Th nucleus. 

According to the recent study Q the energy, E is , of 
the abnormally low-lying level in 229 Th is 7.6±0.5 eV. A 
working compound with an energy gap, A, between the 
top of the valence band and the bottom of the conduction 
band such that A > E is was found in where 232 Th 
nuclei were implanted in LiCaAlFg. There, the Th 4+ ions 
substitute for the Ca 2+ ions and lead to the formation of 
the Th:LiCaAlF 6 alloy with A ~ 10 eV. 

Thus there is a unique situation for low-energy nu- 
clear spectroscopy in 229 Th:LiCaAlF6. There, a photon 
is expected to interact directly with the atomic nucleus 
thereby bypassing the interaction with the electron shells 
5|, i.e. the internal conversion and electronic bridge 
6,0]. This opens up interesting perspectives and ap- 
plications: the 7-ray emission of the excited 229 Th nuclei 
in the optical range called the "nuclear light" Q ; a check 
of the dependence of the isomeric level half life on the 
refractive index [!, @] ; the nuclear metrological standard 
of frequency [|[ called the "nuclear clock" [j| ; the Moss- 
bauer effect in the optical range Q; the relative effect 
of the variation of the fine structure constant e 2 and the 
strong interaction parameter m q / Aqcd (lo| : a check of 
the exponentiality of the decay law [llj . Results obtained 
in the present work can be used in practice to study the 
behavior of the nuclear ensemble with the high density of 
inverse population, the cooperative spontaneous emission 
12] and related phenomena. 

The most important characteristic of the 7-ray laser is 



its amplification factor. The amplification coefficient of 
gamma radiation per unit length, y, caused by stimulated 
emission in a medium is given |l|, [2j by 
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Here we work in the system of units h = c = k = 1, 
Xi S = 2n/E is — 163 ±11 nm is the wavelength of the 
isomeric transition, T ra( i is the radiative width of the 
nuclear transition from the isomeric level to the ground 
state in vacuum, Atu to t is the total width of the transition 
line which accounts for all the types of the homogeneous 
and inhomogeneous broadening, a is the coefficient of the 
electronic conversion, n gr ^ is the density of nuclei in 
the ground (isomeric) state, g — g gr / 9is where g gr (i S ) — 
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gr(is) 



1, and k is the linear attenuation coefficient 



for photons with the wavelength \i S in the medium in all 
processes excluding the resonant transition gr — > is (n ~ 
1 cm^ 1 [13] and it can be reduced). The amplification of 
the 7-radiation occurs if \ > 0. 

The radiating width T ra( i in Eq. ([1} requires the evalu- 
ation of the nuclear matrix element of the Ml transition 
3/2+(7.6 cV)^ 5/2+(0.0). One can calculate this matrix 
element using the Alaga rules and experimental data [lil - 
[3 for the Ml transition 9/2+ (97. 13 keV)->- 7/2+ (71. 82 
keV) Recalculation of the experimental data gives 

respectively the following values for the reduced proba- 
bility of the required transition in terms of the Weisskopf 
units B W . U .(M1) = 0.038, 0.024, and 0.014. The Cori- 
olis interaction between rotational bands enhances the 
transition probability by a factor of 1.2-1.3 [l7j . So for 
the preliminary estimation one can use the "enhanced" 
average value Bw.u.{Ml) ~ 0.032. Correspondingly the 
value of the radiative width is T ra d « 3 x 10 -19 eV. 

The internal conversion of the 3/2 + (7.6 eV) level is for- 
bidden and a — in the 229 Th:LiCaAlF6 crystal because 
A > E is . The probability of the electron bridge [l8| is 
negligibly small in the case A — Ei S > 1 eV [B| . 

In the 229 Th:LiCaAlF 6 crystal the line width Acj tot of 
the transition is determined mainly by the interaction of 
the nuclear magnetic moments [J]. The magnetic mo- 



2 



ment of the ground state fi gr = 0.45/^/v 14l ]. where 
is the nuclear magneton. The magnetic moment of the 
isomeric level /j,i S was calculated in [17j on the basis of 
the available experimental data: a is ~ —0.08u N . With 
these values for u gr and a is the value for the broadening 
of spectral line of the isomeric transition can be reduced 
to 1 kHz Q and approximated by Aw tot < 7 x 10 -13 eV. 

Now we consider the two-step creation of the inverse 
population in the system of the 229 Th nuclei. The first 
step is obtaining a considerable number of the isomers 
229m Th(3/2+, 7.6 eV) by photoexcitation of the 229 Th 
nuclei in the transparent target by the vacuum ultraviolet 
(VUV) laser photons @, H| ■ The equations for n is and 
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as functions of irradiation time are given by 



dn is /dt = <7Lpn gr - A ls n is - gaipn is 



-dng r /dt. (2) 



Here a is the resonance cross section of nuclear photo ex- 
citation by laser radiation with the wavelength \l = \ s 
and the line width Aujl- cr = (Aj s /2ir)T ra d/(gAuJL), f is 
the flux density of laser photons, A is is the decay constant 
of the isomeric level. In our case A is — V ra d — l n (2)/T^ 2 , 
where Tl% « 25 min is the level half-life (we take 1 for the 
refractive index (details see in @, @)). The term gaipn is 
describes the 7 emission of the isomeric nuclei stimulated 
by laser photons. Below we will demonstrate that both 
emission and absorption occur without recoil. Therefore 
the laser radiation is in resonance with the direct and 
return nuclear transitions 5/2+(0.0) +± 3/2+(7.6 eV). 

The solutions of Eq. ((2]) at t — » 00 with the initial con- 
ditions n gr (0) = rig and n.; s (0) = are = noa(p/(Ai S + 
(1 +g)aip), n gr = n (A is + aipg)/{A is + (1 +g)crip). We 
see, that n gr /rii S = Ai S /aip + g > 1, because g > 1. So, 
it is impossible to obtain real inverse population, i.e. to 
fulfill the condition rii s — n gr > by this method. The 
condition m s — n gr /g > in Eq. (fTJ) is not fulfilled too. 

On the other hand, by laser radiation pumping one can 
excite a considerable number of nuclei on the isomeric 
level. As an initial density of the 229 Th nuclei in the 
Th:LiCaAlF 6 crystal we take n = 10 18 cm" 3 0], which 
amounts to 10~ 4 of the LiCaAlFg crystal density. Such 
a small number of the Th 4+ ions in the LiCaAlFe crys- 
tal is not expected to change considerably the electronic 
band structure of the LiCaAlF 6 crystal. Let us consider 
a crystal of 229 Th:LiCaAlFg which has a cylindrical form 
satisfying the condition L 3> D, where D = 0.01 cm is 
the specimen diameter and L is its length. (The Fresnel 
number F = nD 2 / (ALX is ) = 10 if L = 0.5 cm, and F = 1 
if L = 5 cm.) Then for a VUV laser with the average 
power P — 30 mW estimations show that rii s /no ~ 1/4 
during a period of radiation ~ T^%- 

The effective inverse population in the irradiated 
and cooled sample is facilitated by the magnetic field. 
The Zeeman splitting of the nuclear levels is shown in 
Fig. [T] The spin-lattice relaxation time T\ required 
to achieve the Boltzmann distribution is unknown in 
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Figure 1: The transitions between the Zeeman sublevels of 
the 229 Th nucleus in magnetic field. The gray ovals are pro- 
portional to partial level populations. 



229 Th:LiCaAlF6 and should be determined experimen- 
tally. Let us make tentative assumption that Ti meet 
the condition T\ < T^ 2 . In that case populations of the 
sublevels are described by exp(— /x gr ( is )H/T). The ratio 
of the energy differences between the Zeeman sublevels 
of the ground and isomeric states is approximated by 
\ n gr/^is\ ~ 6. Therefore the population of the ground 
state sublevels falls down much faster, than the popu- 
lation of the isomeric sublevels. The magnetic interac- 
tion energy /xH is 1.4 x 10~ 6 eV for the ground state 
and 2.5 x 10~ 7 eV for the isomeric state in the mag- 
netic field H = 100 T. At the temperature T = 0.01 
°K all sublevels of the isomeric state are populated: 
»|3/2+,-a/2> = 9.5 x 10 16 cm" 3 , ^3/2+ ,-1/2) =7.1 x 10 16 
cm" 3 , ni 3 /2+,i/2) = 5.3 x 10 16 cm" 3 , and n| 3 / 2 +,3/2) = 
4.0 x 10 cm -3 . For the ground state calculation gives 
«|5/2+,5/2) = 6.0 x 10 17 cm- 3 , ni 5 / 2 +,3/2> = 1-1 x 10 17 
cm -3 , and n| 5 / 2 +.i/2> = 2.1 x 10 16 cm~ 3 . The popula- 
tion of the other sublevels is negligibly small. 

Thus, the upper Zeeman sublevels of the ground state 
arc unoccupied. Transitions allowed by the selection 
rules for the Ml multipole (they are shown by red ("R") 
arrows in Fig[T]) are possible to these states from the sub- 
levels of the isomeric states. Therefore, in the system of 
nuclear isomers an inverse population is effectively real- 
ized. It is worth noting that there is no resonant ab- 
sorption for the "red" photons because their energies are 
too small to excite the 229 Th nuclei from the populated 
states |5/2+5/2,3/2,l/2}. 

The relative probability p of all these transitions taking 
into account the population of the states and the corre- 
sponding Clebsch-Gordan coefficients is 0.63. This value 
is then used for the calculation of the amplification fac- 
tor, Eq.[TJ where we put n gr = and prii S instead of rii S . 
The result is x — 3 cm -1 . 

The transitions shown by three left green ( "G" ) arrows 
in Fig. Q] also lead to the enhancement of the radiation, 
but their contribution to \ is very small (~ 1%). 

The resonant excitation of the isomeric level from the 
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weakly populated state |5/2 + , 1/2} is also possible (corre- 
sponding transitions are labeled by three right green ar- 
rows in Fig. [TJ) . Such an excitation is a result of the Moss- 
bauer effect in the optical range. Emission of the 7-ray 
photons by the 229m Th isomers and the resonant absorp- 
tion of these photons by the 229 Th nuclei in a solid should 
occur without recoil. Indeed, the energy lost Er due to 
the recoil is negligibly small: E R = uj 2 /2M ps 1.4 x 10" 10 
eV (here M is the 229 Th nucleus mass). The probability 
of the Mossbauer effect at low temperature estimated 
with the Debye-Waller factor is / ~ exp(— 3Er/20e,) 
[l9j . where 9d is the Debye temperature. It is evident 
that in our case / = 1 because Er/8d <C 1. The resonant 
absorption of the gamma radiation is more pronounced 
for the remaining transitions shown in Fig. Q] by three 
left and tree right blue ("B") arrows. 

The sign of fii S (it was predicted theoretically) is not 
critical for the amplification effect. The positive value 
of |Lij s would simply half the general amplification factor. 
In such a case the radiation amplification is due to the 
three transitions to the upper sublevels of the ground 
state from the upper Zeeman isomeric sublevels. 

Let us consider the second scheme of 7-ray laser which 
is based on the electric quadrupole splitting. The 229 Th 
nucleus has appreciable quadrupole moment. The energy 
of the electric quadrupole interaction Eqj is proportional 
to the product of the nuclear quadrupole moment and 
the electric field gradient (EFG), cf> zz . The ground state 
spectroscopic quadrupole moment Q gr is 3.15 eb [lij ]. 
If one makes the standard assumption that the intrinsic 
quadrupole moment Q2 = 8.816 eb [3] remains the same 
for the rotational bands K = 5/2 and K — 3/2 [13] then 
Q zs w 1.8 eb. In the 229 Th:LiCaAlF 6 crystal the leading 
contribution to EFG at the Th 4+ ion site comes from F" 
ions, which compensate the extra charge 2 + . These ions 
are located in interstitial sites in the vicinity of Th 4+ 
2l|,l22j. An estimation gives 4> zz ~ — 10 18 V cm" 2 at the 
Th 4 + site in 229 Th:LiCaAlF 6 . 

To obtain EFG at the 229 Th nucleus one has to take 
into account the effect of antishielding. The Sternheimer 
antishielding factor, 7oo , for Th 4+ is ~100-200. This 
value is extracted from the data approximation [23j for 
7oo m ion series Fr+, Ac 2 +, Ra 3 +, which as Th 4+ have 
the close electron shell of Rn. Thus, EFG at the 229 Th 
nucleus (1 — 7oo)0zz can exceed — 10 20 V cm -2 . Cor- 
respondingly a typical energy of the electric quadrupole 
interaction Eqj can amount to — 10" 4 eV. 

It is also worth noting that according to each 
Th 4+ ion occupies in 229 Th:LiCaAlFg a single lattice site. 
Therefore, the electric quadrupole interaction results in 
constant shift of the transition energy, which does not 
lead to a line broadening. 

If EFG at the 229 Th nucleus in 229 Th:LiCaAlF 6 indeed 
reaches the abovementioned extreme value, then the nu- 
clear splitting scheme will be as shown in FigJ5] (a). The 
sublevels energies are given [2(| by E m — eQ gr u s \(l — 
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Figure 2: The transitions between sublevels of the Th 
nucleus in the case of quadrupole splitting: (a) the case 
Ti < T{ S I2 ; (b) the case T{ a j2 < T x (see text for details). 



loo)4>zz (3m 2 - Jg r (i s ) (J gr (is) + 1))/ (4 Jgr(is) {^Jgr{is) ~ !))• 

At T = 0.01 °K it is mainly the sublevels |5/2+,±l/2) 
(n = 7.4 x 10 17 cm" 3 ) and |3/2+, ±1/2) (n = 2.6 x 10 17 
cm" 3 ) which are populated. The amplifying coefficient 
for the transition |3/2+,±l/2) -)• |5/2+,±3/2) marked 
by the red ("R") arrow in Fig[2ta) is x — 2 cm" 1 . The 
blue ( "B" ) arrows indicate the transitions with the reso- 
nant absorption of the gamma radiation. Notice that at 
large values of EFG the effective inverse population and 
the amplification condition will hold up to the tempera- 
ture T = 0.1 °K. 

Now we consider the third scheme of 7-ray laser when 
TJ% < Ti- In that case the inverse population can be 
achieved by the following way. At first the crystal is hold 
at a low temperature for the time T\ required for the 
Boltzmann distribution population of the sublevels of the 
ground state of the 229 Th nucleus to occur. Thereafter 
the laser radiation provides population of the isomeric 
state. The degree of the crystal heating is governed by 
the absorption coefficient k. Standard purification meth- 
ods can reduce harmful impurities to part-per-billion lev- 
els [j], which ensures the reduction of k to 0.01 cm" 1 . 
An estimation shows that an effective laser heat remov- 
ing even on the scale of few T™ 2 will not increase the 
crystal temperature above T = 0.1 °K. When we stop 
the laser radiation, the cooling down to T = 0.01 °K will 
occur practically instantly. Thus, there will be no re- 
distribution of the nuclear sublevels populations due to 
laser radiation. The upper sublevels of the 229 Th nuclear 
ground state will remain barely populated which will lead 
to final effective inverse population. 

The total cross-section for the isomeric excitations by 
the laser radiation from the |5/2 + ,±l/2) state is vir- 
tually coincides with a in Eq. Furthermore, 90% 
of excited nuclei are in the |3/2 + ,±l/2) state and 10% 
in the |3/2 + ,±3/2) state. The inverse level population 
and the amplification will take place at the "red" transi- 
tions |3/2+,±l/2) |5/2+,±3/2) and |3/2+,±3/2) ->• 
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1 5/2+ ±5/2, ±3/2) in FigHT.). The overall amplification 
factor will be x — 3 cm -1 . 

The duration of the 7-ray laser emission, r, can be 
estimated by the formula r ~ T*J 2 (L/D) 2 exp(— \L) 

X^ 1 ln(ATj S /2), where 



in the case r -C T|% and L < 
N is = n is (nD 2 /4)L. For L = 5 cm r 



10 2 



The 



emission will be a sequence of pulses with the repeti- 
tion frequency f rep — Qi S {D/L) 2 ~ 10 5 -10 4 s _1 , where 
Qi S — Ni S m (2)/T^%. The mean power of the 7-ray laser 
will be ~ 10~ 7 W. For L m 1 cm it is necessary to use sur- 
rounding mirrors to ensure sufficient amplification for the 
light, which makes some passes through the gain medium. 
In this case L in the formula for t should be replaced by 
a total distance which photons pass inside the sample. 

In the following we consider the relaxation time prob- 
lem in detail. According to theoretical models T\ can 
range from few minutes up to millions of years in in- 
sulators at cryogenic temperatures [2qT |. On the other 
hand, there exists a process which can substantially re- 
duce T\. The Ml transitions between the sublevels can 
be accelerated of by covering the sample with a non- 
superconducting metal: Au, Cu, Ag. (Superconductors 
have the energy gap, which can exceed the level split- 
ting AE. It forbids the exchange of the energy AE 
between the nuclei and electrons.) The mechanism of 
the phenomenon is the inelastic scattering of the con- 
duction electrons on the nucleus. Analogous process, i.e. 
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Th reaction, ensures fast decay of the 



6 eV. 
~ 0.5 



Such photon 
h/m* and has the finite range 
0.1 cm. Since r 7 3> D + S, the 



229m Th isomers inside a metal Q_ Let A_E = 10" 
It corresponds to the photon wavelength A = 2n/q 
cm, where the transferred moment q ~ AE^Jm e /2Ep 
(Ep = 5.5 eV is the Fermi energy, m e is the electron 
mass). At low temperature such photon is absorbed 
by a copper layer with depth S < 10 -4 cm. The ex- 
change by the virtual photon between the nuclei inside 
the crystal and the conduction electron in the metal 
cover take place because the effective "mass" of the pho- 
ton to* = \/q 2 — AE 2 ~ q is small 
exists the time At 
r-y ~ cAt — 1/m* 
virtual photon emitted by the nucleus reaches the metal 
cover and interacts with the conduction electrons. This 
process can else be interpreted as internal conversion on 
the conduction electrons. It ensures the nuclear spin re- 
laxation in 229 Th:LiCaAlF 6 . 

The nuclear spin relaxation time is given by \jT\ 
n e {AE / Ep)a e ^DVF-, where n 
conduction electron density, a e ~ 10 - "" cm 
cross-section of the inelastic scattering of the Fermi en- 
ergy electron on the 229 Th ground state sublevels 24|, 
vf = \/2Ep /m e . The factor AE/Ep determines the 
number of electrons near the Fermi surface which can 
absorb or emit the quantum with the energy AE at low 
temperature. Nuclei are placed inside the dielectric sam- 



6 x 10 22 cm 3 is the 



pie with the diameter D. We take into account this fact 
by the factor £d = a e (D)/a e ~ 0.1 (£d can be esti- 
mated through the radial integrals for the electron wave 
functions). Thus, the time required to achieve the Boltz- 
man distribution population of the sublevels of the 229 Th 
ground state is < 50 d. It makes the third scheme of the 
7-ray laser real. It is also worth noting that the metal- 
lic environment of the sample will not affect the 7.6 eV 
transition because A^ s <C D. 

And we close our consideration with a remark on the 
VUV lasers. The energy of the isomeric level is known 
roughly. Therefore we can not tell now, what VUV laser 
will be used for the pumping of the 229m Th isomers. It 
is possible that we can use one of the available lasers 
(see in [IBj]). If not, it will be necessary to develop a 
special laser with the corresponding wavelength or use 
a free electron laser, which has a good tunability. For 
irradiation of the sample we need density of the photon 
flux ip ~ 10 20 cm. _2 s _1 . Such (p can be reached relatively 
easily by focusing of the radiation of middle power laser. 

The author thanks Drs. A. V. Bibikov, I. V. Bodrenko, 
N. N. Delyagin, and A. V. Nikolaev for useful discussions. 
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